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ABSTRACT 


Stream-wise  length  scales  of  coherent  flow  structures  were  examined  to  test  the 
hypothesized  2  to  6  times  the  depth  scale  in  relatively  deeper  (-1  m )  and  faster  (-1.5  m/s) 
flows  in  a  braided  reach  of  the  Kootenai  River,  Idaho.  Velocities  were  measured  using  a 
custom  Acoustic  Doppler  Current  Profiler,  an  array  of  six  electromagnetic  current  meters, 
and  an  Acoustic  Doppler  Velocimeter  in  a  variety  of  river  depths  (0.6-1. 7  m)  and 
velocities  (0.3-1. 6  m/s).  Energetic  (50%  of  total  energy),  coherent  (along  the  array),  low- 
frequency  (<  0.05  Hz)  motions  were  found  for  all  deployment  locations.  Coherent  times 
and  lengths  were  5  s  and  10  m  in  the  ~1.5  m/s  flows  and  45  s  and  22  m  in  the  ~0.5  m/s 
flows.  Multi-resolution  decomposition  provided  coarse  low-frequency  limit  for  the 
coherent  motions  and  suggests  the  temporal  scales  range  from  10  to  1000  s.  Length 
scales  of  the  low-frequency  motions  determined  by  frequency-wavenumber  spectra 
indicate  that  the  motions  are  longer  than  hypothesized.  The  coherent  times  and  lengths 
were  consistently  less  than  the  computed  time  and  length  scales,  suggesting  the  energetic 
low-frequency  motions  evolve  as  they  propagate  downstream. 
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I.  INTRODUCTION 


Flow  in  a  gravel-bed  river  appears  to  be  random  and  chaotic,  but  through 
laboratory  and  field  experiments,  it  has  been  shown  that  it  is  possible  to  classify  the 
flow  into  organized  motions  (Cantwell  1981;  Roy  et  al.  2004),  which  possess  spatial 
and  temporal  coherence  (Adrian  2007).  An  understanding  of  river  flow  characteristics 
is  important  for  describing  river  morphology,  bed-load  transport  (Robert  et  al.  1996), 
and  mixing.  River  flow  is  generally  divided  into  three  vertical  layers:  1)  the  near- wall 
bottom  boundary  consisting  of  less  than  10%  of  the  flow  depth;  2)  the  outer  flow  region 
containing  the  largest  percentage  of  flow;  and  3)  the  near-surface  region  (Tamburrino 
and  Gulliver  1999).  The  two  main  mechanisms  believed  responsible  for  generating 
large-scale  coherent  flow  structures  are  turbulence,  as  a  result  of  bed  topography  in  the 
near- wall  region  interacting  with  the  outer  flow  region  (Hardy  et  al.  2010),  and 
coalescence  of  vortices  generated  from  prominent  river  structures  (Cantwell  1981). 
Near-bed  turbulence  can  produce  coherent  flow  structures  during  a  bursting  (upwelling) 
event,  where  low  speed  fluid  is  ejected  into  the  outer  flow  region  that  persist  until  a 
sweeping  (downwelling)  event  returns  high  speed  flow  energy  to  the  near-wall 
boundary  layer  (Hardy  et  al.  2009).  In  the  wake  of  river  features,  such  as  large  pebble 
clusters,  riffle  pools,  or  protruding  structures,  vertical  and  horizontal  eddies  are 
produced  at  a  constant  rate  and  combine  to  produce  large-scale  coherent  flow  structures 
(Buffin-Belanger  and  Roy  1998;  Cantwell  1981).  In  both  cases,  the  structures  have 
been  described  as  recurrent  patterns  of  alternating  high  and  low  speed  wedges  (Buffin- 
Belanger  et  al.  2000),  repetitive  quasi-cyclic  large-scale  turbulent  motions,  or  simply, 
rolling  vortices  that  move  downstream  at  flow  velocity  (Shvidchenko  and  Pender 
2001). 

Coherent  flow  structure  size  is  believed  to  be  controlled  by  flow  depth,  and 
experiments  using  flow  visualization  techniques  and  field  measurement  analysis  usually 
provide  size  scales  as  a  function  of  flow  depth  (e.g.,  Shvidchenko  and  Pender  2001; 
Roy  et  al.  2004).  Since  temporal  coherence  is  a  prerequisite  for  identifying  a  coherent 
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flow  structure  (Adrian  2007),  the  length  scale  is  commonly  estimated  using  Taylor’s 
frozen  turbulence  hypothesis  in  which 


L  =  UT,  (1) 

where  U  is  the  mean  flow  velocity  and  T  is  the  time-scale,  or  average  period,  of  the 
motion  (Soulsby,  1980).  Zaman  and  Hussain  (1981)  concluded  Taylor’s  hypothesis  is 
acceptable  for  application  to  large-scale  coherent  flow  structures.  Most  studies 
consistently  find  that  the  vertical  scale  of  coherent  flow  structures  is  proportional  to 
flow  depth,  but  the  horizontal  length  scales,  at  least  in  the  stream-wise  direction,  have 
not  been  precisely  defined.  Most  studies  provide  a  mean  lateral,  or  cross-stream,  width 
scale  of  1  to  2 h,  where  h  is  the  flow  depth,  and  a  mean  longitudinal,  or  stream-wise, 
length  scale  of  2  to  6 h  (Marquis  and  Roy  2006).  Shvidchenko  and  Pender  (2001)  using 
flow  visualization,  quadrant  analysis,  and  energy  estimates  reported  a  mean  stream- 
wise  length  scale  of  4  to  5 h  for  coherent  flow  structures  in  a  laboratory  flume  study. 
Roy  et  al.  (2004),  amongst  others,  used  several  techniques  described  by  Nakagawa  and 
Nezu  (1981)  on  field  measurements  acquired  in  rivers  with  similar  velocities  (<0.5  m/s) 
and  depths  (<0.5  m)  and  found  a  3  to  5 h  length  scale.  Dinehart  (1999)  conducted  one 
of  the  few  known  coherent  flow  structure  studies  in  a  gravel-bed  river  with  a  relatively 
fast  (2  m/s)  and  deep  flow  (1.5  m)  using  14.3  m  instrument  separation.  He  found 
coherent  times  of  10  to  30  s  with  lengths  of  coherent  flow  structures,  termed  correlated 
stream- wise  eddies,  ranging  from  20  to  30  m  that  correspond  to  a  13  to  20/?  scale. 
These  longer  length  scales  for  coherent  flow  structures  in  a  deeper  and  faster  flow  and 
the  range  of  postulated  length  scales  from  others,  points  out  the  need  for  a  detailed 
study  in  which  measurements  are  collected  in  a  range  of  depths  and  flow  velocities. 

Two  common  techniques  to  identify  coherent  flow  structures  are  deploying 

either  synchronous  lagged  instrument  arrays  (e.g.,  Roy  et  al.  2004)  or  separated 

individual  instruments  placed  some  distance  apart  (e.g.,  Dinehart  1999).  Both 

techniques  measure  stream-wise  velocities  and  various  methods  are  applied  to  estimate 

the  coherence  between  sensors  associated  with  the  passage  of  a  flow  structure.  A 

synchronous  lagged  instrument  array  is  required  to  describe  coherent  flow  structures, 

but  this  requires  cabling  and  limits  the  array  length.  It  is  also  limited  by  the  number  of 
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expensive  sensors  in  the  array,  sensor  placement,  and  the  requirement  of  a  frame 
assembly  to  mount  the  sensors,  cables,  data  logger,  and  power  source  in  the  river  at 
fixed  locations  while  minimizing  flow  disturbance.  These  limitations  have  constrained 
previous  studies  to  laboratory  or  shallow  water  field  experiments  that  allow  for  shorter 
arrays. 

The  objective  of  this  study  is  to  measure  stream- wise  length  scales  of  large- 
scale  coherent  flow  structures  in  deeper  and  faster  flow  conditions  of  a  gravel-bed  river. 
This  study  was  conducted  as  part  of  a  research  initiative  sponsored  by  the  Office  of 
Naval  Research  (ONR)  to  investigate  the  sensitivity  of  Delft3D  river  models  to 
bathymetric  variability.  The  purpose  is  to  increase  understanding  of  river  flow 
characteristics  to  improve  model  output.  Instruments  were  deployed  in  various  flow 
depths  and  velocities  to  investigate  the  effect  of  these  variables  on  the  purported  length 
scales  of  coherent  flow  structures.  This  study  considers  coherent  motions  as  periodic 
structures  resembling  a  wave  packet,  similar  to  Adrian’s  (2007)  eddy  packet 
description,  that  propagate  in  the  stream-wise  direction  and  applies  Fourier  based 
methods  for  isolating  motions  within  the  flow.  The  capability  of  an  Acoustic  Doppler 
Current  Profiler  (ADCP)  oriented  with  a  horizontal  beam  in  the  stream-wise  direction  is 
evaluated  to  non-intrusively  measure  flow  velocities  at  large  distances.  An  Iterative 
Maximum  Likelihood  Estimator  (IMLE)  method  is  used  that  can  resolve  length  scales 
up  to  two  or  four  times  the  maximum  array  length,  which  has  not  been  previously  used 
in  riverine  environments.  The  results  will  show  that  the  resolution  of  data  analysis 
techniques  is  an  important  factor  in  providing  accurate  delineation  of  the  motions,  and 
thus  length  scales,  present  in  the  flow. 
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II.  EXPERIMENT  DESIGN 


A.  EXPERIMENT  SITE  DESCRIPTION 

River  velocities  were  acquired  by  three  different  sensors  as  part  of  an  8-day 
riverine  field  experiment  conducted  in  August  2010  on  a  braided  reach  of  the  Kootenai 
River  near  Bonners  Ferry,  Idaho.  The  Kootenai  River  flows  south  from  its  headwaters 
in  the  Rocky,  Salish,  and  Purcell  Mountain  ranges  through  northwestern  Montana,  turns 
north  and  flows  through  the  Kootenai  valley  of  northern  Idaho,  and  reaches  its  terminus 
at  Kootenay  Lake  in  British  Columbia,  Canada.  The  discharge  of  the  Kootenai  River  is 
primarily  controlled  by  the  Libby  Dam  in  Montana  and  was  a  constant  196  m  /s  during 
the  experiment  (U.S.  Geological  Survey  2010).  The  experiment  site  is  located  on  two 
relatively  straight  120  m  wide  sections  of  the  main  river  channel  and  a  narrower  20  m 
wide  straight  side  channel  (Figure  1).  Bed  substrate  ranged  from  coarse  gravel  to 
cobbles  with  particle-size  diameters  between  20.8  and  78.2  mm  ( Di6  and  Dx4),  with  a 
median  particle  size  (D^o)  of  39.9  mm,  where  Dx  represents  the  diameter  of  the  xth 
percentile  from  the  bed  particle  distribution  (Fosness  and  Williams  2009).  The  river 
bank  is  comprised  of  rough  gravel  and  is  approximately  3  m  wide  followed  by  a  steep  6 
m  elevation  rise  to  the  floodplain. 

The  site  is  divided  into  four  distinct  zones  characterized  by  river  feature  and 
flow  pattern  (Figure  2).  Zone  1  encompasses  a  section  of  the  main  channel  with  a 
thalweg  depth  >1  m  and  mean  flow  velocity  >1.6  m/s.  Zone  2  has  similar  flow 
characteristics  as  Zone  1  and  is  located  2  km  upstream.  Zone  3  includes  the  20  m  wide 
side  channel  with  a  thalweg  depth  of  1.6  m  and  mean  flow  velocity  of  0.5  m/s.  The 
channel  narrows  and  shallows  to  0.7  m  at  a  50  in  long  riffle  section  causing  the  flow  to 
increase  to  1.4  m/s.  Zone  4  is  located  downstream  of  a  man-made  concrete  piling  100 
m  upstream  of  Zone  3  with  a  mean  depth  and  velocity  of  1  in  and  0.4  m/s.  Downstream 
of  the  structure,  the  flow  splits  at  a  gravel  bar  and  either  continues  through  to  the  side 
channel  or  re-joins  the  main  channel,  where  the  confluence  affects  the  flow  pattern  of 


5 


Zone  1  near  a  large  gravel  bar.  Zones  1  and  2  are  representative  of  typical  river  flows 
and  Zones  3  and  4  contain  areas  of  visible  vortex  shedding. 


B.  INSTRUMENTATION 

A  synchronous  lagged  array  scheme  was  implemented  to  compute  the  expected 
stream- wise  length  scales  of  coherent  flow.  Stream- wise  velocities  were  measured 
using  a  non-equidistant  lagged  array  of  six  Valeport  electromagnetic  current  meters 
(ECM),  a  custom  2  MHz  Nortek  ADCP,  and  a  Nortek  Acoustic  Doppler  Velocimeter 
(ADV)  mounted  on  a  GPS-equipped  portable  aluminum  frame  (Figure  3).  The  frame 
was  designed  to  support  the  4  m  long  ECM  array,  oriented  in  the  stream-wise  direction, 
and  to  minimize  flow  disturbances.  The  frame  was  constructed  of  5  cm  diameter  pipes 
and  resembled  a  2  m  square  scaffold.  The  bracing  and  legs  were  adjustable  to  locate 
and  level  the  frame  in  a  variety  of  river  bottom  contours.  The  spacing  between  each  6.5 
mm  diameter  ECM  was  set  to  maximize  the  resolution  of  coherent  motions  from  0.5  m 
up  to  8  m  in  length  using  an  IMLE  method  described  later.  The  ADCP  was  fitted  with 
a  custom  head  (Appendix  A)  that  uses  one  horizontal,  1.7  degree  beam  width  to 
measure  along  beam  velocities  in  the  stream-wise  direction  at  equal  35  cm  bins,  up  to 
10  m  in  length,  with  a  sampling  frequency  of  1  Hz.  The  stream-wise  oriented  beam 
conveniently  serves  as  a  high-resolution,  equidistant  lagged  instrument  array  without 
the  limitations  of  typical  synchronous  arrays,  such  as  cabling,  a  large  frame,  and  a 
separate  data  logger.  To  date,  this  is  the  first  experiment  to  use  an  ADCP  to  compute 
scales  of  coherent  flow  structures.  Since  ADCPs  are  inherently  noisy  instruments,  and 
the  perfonnance  of  ADVs  and  ECMs  in  measuring  velocities  in  turbulent  flows  (~1.5 
m/s)  has  been  shown  suitable  for  describing  flow  characteristics  (MacVicar  et  al.  2007), 
ADCP  velocities  are  compared  to  the  ECM  and  ADV  (Appendix  A). 

Stream-wise  velocities  were  measured  for  a  minimum  of  1.5  hours  at  each  site. 
The  ADCP  and  ADV  were  mounted  approximately  1.5  m  apart  in  the  across-stream 
direction  on  the  front  of  the  aluminum  frame  and  all  sensor  heads  were  placed  at 
approximately  mid-depth  (Figure  3).  The  frame  was  deployed  in  several  locations  of 
the  river  (Figure  1  and  Table  1)  that  included  the  large,  fast  flow  main  channel,  the 
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smaller,  slower  flow  side  channel,  downstream  of  man-made  structures,  and  upstream 
and  downstream  of  riffles.  The  scheme  for  instrument  deployment  was  systematic  for 
Zones  1,  2  and  3,  starting  near  the  river  bank  at  approximately  0.5  m  depth  and 
progressively  moving  to  faster  and  deeper  flows  farther  out  towards  mid-channel.  In 
Zones  1  and  2,  however,  the  strong  stream-wise  current  (-1.6  m/s)  near  mid-channel 
almost  caused  the  instrument  frame  to  topple  and  prevented  velocity  measurements  for 
depths  greater  than  1  m  for  these  areas  of  the  river.  The  slower  velocities  in  Zones  3 
and  4  did  not  limit  frame  deployments  in  depths  up  to  1.5  m. 

C.  DATA  VALIDATION 

Initial  data  verification  of  the  first  deployment  showed  that  the  ECM  array 
measurements  were  affected  by  electromagnetic  interference  from  overhead  power 
lines  and  these  datasets  were  removed  from  analysis.  The  ECM  array  was  then 
deployed  farther  upstream  away  from  interference  in  Zone  2  for  two  days  to  provide  a 
benchmark  for  ADCP  validation  (Appendix  A).  The  larger  number  of  ADCP  datasets 
and  longer  beam  lengths,  compared  to  the  ECM  array,  made  these  measurements  ideal 
for  use  in  all  lagged  array  analysis  methods  in  this  study.  Several  other  datasets  were 
not  included  in  the  analysis  due  to  instrument  orientation  where  either  the  effective 
beam  length  of  the  ADCP  was  too  short  for  analysis,  or  river  features,  such  as  riffles, 
prevented  accurate  coherent  flow  structure  identification.  In  all,  measurements  from  10 
of  the  initial  19  deployment  locations  were  used  for  analysis.  A  three  standard 
deviation  filter  to  remove  outliers  was  applied  to  the  data  (Emery  and  Thomson  2001). 
For  all  ADCP  datasets,  the  mean  velocity  and  energy  spectra  per  bin  were  plotted  to 
identify  the  maximum  bin  where  the  mean  velocity  and  spectral  curve  deviated  more 
than  10%  from  the  prior  bin’s  measurement. 
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III.  DATA  ANALYSIS  METHODS  AND  RESULTS 


A.  SPECTRAL  ANALYSIS 


Spectra  were  calculated  by  dividing  the  1.5-hour  velocity  time  series  into  20 
Hanning  windowed  subseries  of  256  5  with  50%  overlap,  resulting  in  a  frequency 
resolution  of  0.0039  Hz  and  106  degrees  of  freedom.  The  process  to  detennine  ideal 
record  and  subseries  length  is  explained  in  Appendix  B.  Longer  (>3  hours)  ADCP 
records  at  fixed  locations  were  divided  into  several  1.5-hour  record  lengths  providing 
15  datasets  for  analysis.  The  energy  spectra  are  qualitatively  separated  into  frequency 
bands  and  used  to  compute  the  respective  scales.  The  average  period  for  the  motions  in 
each  frequency  band  is  computed  by 


[%(/¥/ 

J  J 1 _ 


(2) 


where  E  is  the  energy  density  and/ is  the  frequency. 


Spectral  estimates  of  representative  deployment  locations  for  each  zone  are 
shown  in  Figure  4.  The  spectra  consistently  show  energy  in  the  lower  frequencies  with 
a  well-defined  energy  cascade  occurring  between  -0.02  to  0.1  Hz.  A  noise  floor  is 
evident  and  consistent  in  all  spectra  with  an  amplitude  between  0.01  to  0.02  (m/s)  -s. 
Zone  3  and  4  spectra  show  distinct  and  significant  peaks  at  0.01  Hz  indicating  the 
presence  of  a  strongly  periodic  motion  and  is  attributed  to  the  formation  of  vortices  that 
are  shed  from  the  riffle  (Zone  3)  and  the  man-made  concrete  piling  (Zone  4). 
Horizontal  eddies  were  visually  observed  downstream  from  these  features  and  along  the 
high-speed  jet  of  fluid  ejected  from  the  riffle.  Application  of  Taylor’s  hypothesis  (Eq. 
1)  to  the  low-frequency  periodic  motions  in  Zones  3  and  4  results  in  a  temporal  and 
length  scale  of  100  s  and  44  to  53  m,  however,  it  is  unclear  whether  the  motion  applies 
to  the  horizontal  eddies,  vertical  flow  structures,  or  their  combination.  An  anomalous 
spectral  peak  (not  shown)  was  apparent  in  a  Zone  2  deployment,  which  had  a  nearly  5 
hour  total  record  length  and  was  divided  into  three  smaller  1.5-hour  record  lengths. 
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This  anomaly  prompted  an  investigation  and  test  into  smaller  record  lengths  (Appendix 
B),  and  the  test  revealed  motions  that  appear  periodic  in  the  spectra  from  record  lengths 
shorter  than  1.5  hours  are  not  significantly  different  than  random  noise.  Since  the 
peaks  in  Zones  3  and  4  are  from  energy  spectra  computed  with  1.5-hour  ADCP  record 
lengths,  the  motions  are  considered  significant.  An  energy  roll-off  towards  zero  is  not 
apparent  in  the  lowest  frequency  bin  in  Zones  1  and  2.  Similar  results  were  obtained  by 
varying  the  subseries  record  lengths  to  increase  low-frequency  resolution  (Figure  5). 
This  suggests  that  longer  record  lengths  (>  1.5  hours)  are  needed  to  increase  low- 
frequency  resolution  while  maintaining  statistical  confidence.  In  many  cases,  the 
lowest  frequency  bin  approaches  zero,  but  the  increased  resolution  masks  any  apparent 
energy  cascade  and  separation  of  motions  is  difficult.  An  accurate  temporal  scale  for 
the  low-frequency  motions  in  Zones  1  and  2  cannot  be  adequately  determined  from  Eq. 
2  without  complete  frequency  limits. 

Normalized  logarithmic  spectra  for  each  ADV  deployment  were  plotted 
together  (Figure  6)  and  linear  regression  used  to  compute  slopes  associated  with  the 
features  visible  in  the  spectra.  The  motions  in  the  low-frequency  band  (<  0.05  Hz)  have 
a  mean  slope  near  zero  suggesting  the  presence  of  large-scale,  low-frequency 
structures.  The  roll-off,  or  cascade,  of  energy  in  the  -0.02  to  0.1  Hz  frequency  band  is 
typical  of  turbulent  motions,  but  the  mean  slope  of  -1  differs  from  the  expected  -5/3 
slope  associated  with  turbulence.  Similar  spectral  slopes  were  noted  in  a  laboratory 
experiment  on  turbulent  river  flows  conducted  by  Carrasco  and  Vionnet  (2004).  This 
suggests  the  presence  of  an  inverse  energy  cascade,  where  energy  is  transferred  to  the 
low-frequency  motions  from  higher  frequency  turbulence.  An  energy  cascade,  or 
spectral  gap,  usually  detennines  the  separation  of  motions  present  in  the  flow  (Carrasco 
and  Vionnet  2004;  Vickers  and  Mahrt  2003)  and  based  on  experiments  by  others  (e.g., 
Sukhodolov  and  Uijttewaal  2010),  the  energy  cascade  represents  turbulence  within  the 
flow  and  separates  the  low-frequency,  energetic  motions  from  the  noise  floor.  The 
frequency  band  of  turbulence  is  qualitatively  determined  and  noted  as  an  upper  limit  for 
the  low-frequency  motions.  The  limit  ranges  between  0.0078  to  0.0625  Hz  and  is 
consistent  for  all  datasets  of  the  same  deployment.  Roy  et  al.  (2004)  noted  a  frequency 
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band  of  0.07  to  0.12  Hz  for  coherent  flow  structures  in  a  relatively  slow,  shallow  flow 
(0.67  m/s  and  0.32  m).  They  used  a  20  min  total  record  length  divided  into  60  5 
subseries  lengths.  Their  choice  of  window  length  is  described  in  Buffin-Belanger  and 
Roy  (2005)  in  which  the  low-frequency  spectra  contained  <3  %  of  the  total  signal 
energy  (variance),  and  therefore,  in  their  study,  the  record  lengths  could  be  truncated  to 
disregard  the  low-frequency  motions. 

Normalized  cumulative  integrated  spectra  were  computed,  which  represents  the 
percentage  contribution  of  the  motions  to  the  total  variance  per  frequency  (Figure  7). 
Disregarding  the  noise  floor  (f  >  0.25  Hz),  the  spectra  show  that  between  40  to  50%  of 
the  variance  is  contained  in  the  low-frequency  (<0.05  Hz)  band  prior  to  the  energy 
cascade.  This  result  indicates  that  the  low-frequency  motions  in  this  study  are 
significant  and  cannot  be  ignored. 

B.  COHERENCY 


Signal  coherency  is  a  Fourier  based  method  to  measure  the  linear  correlation  by 
frequency  of  two  signals  acquired  by  two  sensors  with  spatial  separation.  Two  time 
series  are  transfonned  into  the  frequency  domain  through  Fourier  techniques  and 
compared  to  detennine  if  the  signals  are  coherent  as  a  function  of  frequency  (Emery 
and  Thomson  2001).  The  squared  coherency  function  between  two  sensors  is  defined 
as 


=  cixn+QUf) 

"  K1(/)||s2,2(/)| 


(3) 


where  Sn(f)  and  S22(f)  are  the  auto-spectra  of  the  two  signals,  C“  (/) 


and 


tfc f) 


are  the  co-spectrum  and  quadrature  spectrum  of  the  cross-spectrum,  and 


0<  -  !•  Two  completely  correlated  signals,  without  noise,  will  have  a  y 1 


equal  to  1  while  two  un-correlated  signals  will  have  a  y1  equal  to  0.  A  y2  between  0 
and  1  occurs  if  either  the  signals  are  not  linearly  related  or  noise  is  contained  within  the 
measurements  (Bendant  and  Piersol  2000).  Confidence  limits  for  y2  are  specified  to 
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delineate  a  significance  level  threshold  at  which  point  the  signals  are  no  longer 
statistically  similar.  A  significance  level  of  0.06  is  obtained  for  a  confidence  level  of 
95%  from  a  chi-square  distribution  with  106  degrees  of  freedom  (Emery  and  Thomson 
2001).  Therefore,  two  velocity  signals  are  considered  coherent  if  y2  is  between  0.06 
and  1.  The  maximum  sensor  separation  that  retains  coherency,  as  a  function  of 
frequency,  is  the  coherent  length  ( Lc )  of  the  signal.  Lc  indicates  the  distance  at  which 
the  signal  is  statistically  the  same  or,  in  other  words,  the  distance  that  the  signal  persists 
before  it  can  no  longer  be  considered  the  same  signal.  The  Lc  of  a  signal  is  not  an 
indication  of  its  total  motion  length,  as  a  signal  can  persist  longer  or  shorter  than  its 
length.  The  coherence  time  ( Tc),  which  is  the  persistence  time  of  a  coherent  signal,  is 
calculated  from  Lc  by 

(A 

c  U 

assuming  Taylor’s  hypothesis.  This  technique  is  used  to  confirm  the  presence  of 
coherent  motions,  as  a  function  of  frequency,  and  to  provide  Lc  and  T,  of  the  large  scale 
motions. 

Stream-wise  velocity  coherency  was  computed  at  the  first  ADCP  bin  relative  to 
each  successive  bin  up  to  the  maximum  beam  length  (Figure  8  a-b).  Coherency 
decreases  to  lower  frequencies  with  increasing  lag  suggesting  larger  Lc  for  the  low- 
frequency  motions.  The  maximum  lag  at  which  coherency  is  retained  ( y2>  0.06)  is 
determined  and  plotted  as  a  function  of  frequency  (Figure  9).  The  maximum  distance  is 
long  for  low  frequencies  and  decreases  rapidly  for  higher  frequencies.  A  triangular  plot 
marker  is  used  to  indicate  the  instances  when  the  maximum  lag  at  which  coherency 
occurs  equals  the  maximum  bin  size.  Lc  for  the  low-frequency  motions  are  longer  than 
the  beam  lengths  of  the  representative  flows  measured  in  each  zone  (Figure  9).  In 
Zones  1  and  3,  the  low-frequency  coherent  motions  persist  farther  than  6  to  9  m.  The 
mean  y2  for  the  energetic,  low-frequency  band  motions  in  Zones  1  and  2  and  y2(f  = 
0.01  Hz)  for  Zones  3  and  4  appear  to  linearly  decrease  as  a  function  of  bin  lag  (Figure 
10).  A  linear  best-fit  line  suggests  the  low-frequency  coherent  motions  persist  up  to  10 
m  (Zones  1  and  2)  and  22  to  27  in  (Zones  3  and  4).  Finearly  extracted  Lc,  with 
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corresponding  Tc  and  Lc/h  values  are  listed  in  Table  2.  Low-frequency  energetic 
motions  are  at  least  temporally  coherent,  have  relatively  long  Lc,  and  therefore,  this 
study  considers  these  motions  as  large-scale  coherent  flow  structures. 

C.  FREQUENCY- WAVENUMBER  ANALYSIS 

Next,  the  stream-wise  spatial  scales  for  the  coherent  flow  structure  are  extracted 
from  frequency-wavenumber  (j-ks)  spectral  analysis.  This  method  has  been  used 
extensively  in  nearshore  experiments  to  compute  long,  infragravity  wave  energy  as  a 
function  of  frequency  and  wavenumber  (Pawka  1983;  Oltman-Shay  and  Guza  1987; 
Oltman-Shay  et  al.  1989;  Reniers  et  al.  2006).  These  periodic  waves  have  temporal  and 
spatial  scales  of  20  to  200  s  and  20  to  1000  m  (Oltman-Shay  and  Guza  1987).  The 
wavenumber  (inverse  length),  k,  is  a  vector  comprised  of  an  across-stream  component 
( ka )  and  a  stream-wise  component  (ks).  Lagged  sensor  arrays  are  usually  deployed,  and 
f-ks  spectra  for  each  frequency  are  computed  for  propagating  wave  components  (Capon 
1969).  The  f-ks  spectra  are  detennined  from  the  spatially  lagged  cross-spectrum 
[M(£,f)\  by 

E(f,ks)  =  j*M(f,/)exp(-iA:s  *i)dl ,  (5) 

where  E(f,ks)  is  the  spectral  value  at  frequency / and  stream-wise  wavenumber  ks  and  t 
is  the  spatial  lag.  The  complex  cross-spectrum  is  determined  by  Fourier  decomposition 
according  to 

Mij{f)  =  Ff{v(si,t)}F}{v(sj,t)},  (6) 

where  Ff{  }  is  the  complex  Fourier  coefficient  at  frequency,  f  of  the  velocity 
component,  v,  at  the  stream-wise  locations  Sj  and  Sj.  Two-dimensional  Fourier 
decomposition  to  produce  the  cross-spectrum  matrix  [My(f)\  can  only  be  computed  at 
the  discrete  equal  lag  separations  of  the  sensors  in  the  array  which  yields  low  resolution 
results  for  E(fks)  (Sheremet  et  al.  2005). 

If  the  motion  lengths  of  interest  are  longer  than  the  array  length  or  the 
instrument  array  consists  of  unequally  spaced  lags,  high-resolution  estimation  methods 
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can  be  used  (Oltman-Shay  and  Guza  1987).  An  IMLE  method  provides  a  linear 
estimator  of  E(f,ks)  and  has  the  form 

E{f,K)  =  YLa»V’K)My<J),  (V) 

i  j 

where  (/,  ks )  are  complex  variable  weights  composed  of  sines  and  cosines  (Pawka 

1983).  The  cross-spectra  are  computed  at  each  frequency,  then  IMLE  applies  a  best- fit 
curve  algorithm  to  the  spectra,  and  ks  values  are  extracted  from  the  curve  fit.  This 
curve  fit  is  optimal  when  at  least  two  points  of  the  passing  wave  are  captured  in  the 
cross-spectrum.  Tests  perfonned  by  Oltman-Shay  and  Guza  (1984)  show  the  IMLE 
method  provides  the  most  accurate  estimate  of  the  energy  contained  per  wavenumber 
bandwidth,  and  total  energy  estimates  were  accurate  to  within  1%  of  the  true  spectrum. 
IMLE  extends  the  array  by  extrapolating  the  measurements  and  increases  the 
wavenumber  resolution  up  to  four  times  the  array  length  with  more  statistically 
confident  estimates  at  twice  the  array  length  (Sheremet  et  al.  2005).  The  co-spectra  for 
various  frequencies  of  a  Zone  2  deployment  are  shown  in  Figure  1 1  to  emphasize  the 
IMLE  process.  In  the  low  frequencies  (<0.063  Hz),  less  than  a  half  motion  length  is 
computed  and  therefore,  wavenumbers  cannot  be  confirmed  below  this  frequency.  The 
co-spectra  at  a  higher  frequency  (0.49  Hz)  reveal  that  the  motions  decay  at  longer  lags. 
This  is  a  common  result  of  this  technique  (e.g.,  Oltman-Shay  et  al.  1989,  Figure  5)  and 
is  similar  to  the  results  of  coherency  analysis,  since  both  techniques  use  co-spectra 
computations  in  their  formulations.  The  IMLE  method  is  used  to  compute  the 
estimated/-^,  spectra  for  the  stream-wise  velocities  (Figure  12). 


A  length  scale  is  obtained  from 

l  -  <  SA _  (i 

k)  t,j;  E(f,k,)dfdk,’ 

where  Ay  is  the  wavenumber  limit  detennined  by  two  or  four  times  the  maximum  array 
length  and  Ay  is  the  Nyquist  wavenumber  determined  by 


1 

2M  ’ 
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where  M  is  the  smallest  instrument  lag  separation  in  the  array.  The  wavenumber 
limit,  k  ,  determines  the  resolvable  limit  of fj.  Wavenumbers  computed  by  doubling 

the  array  length  provide  statistically  conservative  estimates,  while  those  computed  with 
four  times  the  array  still  provide  reliable,  but  less  confident,  estimates.  Length  scales 
are  obtained  by  averaging  the  wavenumbers  over  the  frequency  bands  (fi-fi)  identified 
from  energy  spectra  that  separate  the  low-frequency,  coherent  flow  structures  and 
turbulence. 

A  ridge  of  energy  as  a  function  of f-ks  follows  the  advective  velocity  of  the  flow 
according  to  Taylor’s  hypothesis  (Eq.  1),  highlighted  in  Figure  12.  The  f-ks  spectra 
show  concentrations  of  energy  in  low-frequency,  small  wavenumber  (large  motion 
length)  bands  suggesting  the  coherent  flow  structures  have  relatively  long  periods  and 
lengths.  The  spectra  reveal  that  the  wavenumbers  for  the  low-frequency,  coherent  flow 
structures  are  beyond  the  resolvable  range  of  the  array  length.  In  the  spectra  of  a  Zone 
3  deployment,  with  a  beam  length  of  9  m,  a  small  frequency  band  of  the  coherent  flow 
structure  was  available  to  extract  a  length  of  around  7  m  (Table  3).  The  actual  mean 
structure  length  should  be  much  larger  and  could  be  confirmed  if  the  ADCP  beam 
length  was  longer,  but  nevertheless,  this  short  structure  length  is  larger  than  the 
predicted  2  to  6 h  scale.  Mean  motion  lengths,  L,  averaged  over  the  turbulence 
frequency  band,  are  computed  from  Eq.  8  for  each  dataset  (Table  3).  In  some  datasets, 
indicated  by  an  asterisk  in  Table  3,  the  array  length  prevented  the  detennination  of  the 
full  range  of  motion  lengths  for  the  turbulence.  The  actual  lengths  for  the  low- 
frequency,  coherent  flow  structures  are  much  larger  than  anticipated,  highlighting  the 
need  for  a  longer  (»  10  m)  array  length. 

D.  MULTI-RESOLUTION  DECOMPOSITION 

To  address  the  resolution  of  a  lower  frequency  limit  for  the  coherent  flow 
structures,  multi-resolution  decomposition  (MR)  has  been  used  in  analyzing  turbulent 
fluxes  and  provides  coarse  time-scales  of  fluctuations,  or  events  (Vickers  and  Mahrt 
2003).  MR  is  a  type  of  spectral  averaging,  which  computes  sample  variance  associated 
with  2P  averaging  length  time-scales  and  represents  the  simplest  orthogonal 


15 


decomposition  (Vickers  and  Mahrt  2003).  Howell  and  Mahrt  (1997)  describe  the 
spectra  produced  from  MR  as  simultaneous  decomposition  of  a  signal  into  un-weighted 
moving  averaging  windows  of  different  widths.  Both  Fourier  methods  and  MR 
preserve  signal  variance,  but  unlike  Fourier  where  the  spectral  peak  in  the  frequency 
domain  depends  on  periodicity,  the  peak  of  MR  spectra  in  the  time  domain  depends  on 
time-scale  fluctuations  in  the  signal  (Vickers  and  Mahrt  2003).  A  complete  explanation 
of  MR  is  given  in  Howell  and  Mahrt  (1997)  and  a  brief  description  of  the  MATLAB 
software  algorithm  provided  by  Vickers  (2010,  personal  communication)  is  presented 
herein.  Given  a  time  series  of  2P  data  points,  MR  first  partitions  the  record  into  simple 
averages  of  different  segments  of  width  1,  2,  4.  ,.2P  points.  The  segments  are 
sequenced  as  n  =  2P,  2P'  ...1,  where  n  is  the  position  of  the  segment  within  the  series. 
The  averages  are  removed  from  each  consecutive  sub-segment  starting  with  the  largest 
scale  ( 2P  width)  and  ending  with  the  shortest  time-scale.  The  average  for  the  nth 
segment  of  scale  p  is  given  by 

v«0)  =  Ty  £  vr(p),  (TO) 

2-  ;=(n-l)2p+l 

where  vr/p)  is  the  residual  series  after  segment  averages  with  windows  of  width  >  2P 
points  are  removed.  The  process  starts  with  the  largest  segment  (the  entire  time  series) 
and  removes  the  mean  from  all  the  data  points.  The  resulting  time  series  is  the  1st 
residual  series.  The  next  step  moves  to  the  next  smaller  segment,  which  is  half  the 
length  of  the  original  time  series  and  composed  of  1st  residual  series  values  and 
removes  the  segment  mean  from  the  respective  segment  data  points  resulting  in  the  2nd 
residual  series.  The  process  continues  to  remove  the  segment  means  from  each 
segment  as  it  proceeds  through  to  the  shortest  time-scale  associated  with  the  record 
length.  The  average  for  each  nth  segment  is  used  to  calculate  the  MR  spectra,  which 
represents  the  variance  of  the  segment  averages  and  is  given  by 

A(/>+1)  =  AtZ7(/’).  (11) 

^  n=\ 

where  vn(p)  is  the  segmented,  mean  removed  time  series,  Dv(p+1)  is  the  variance  of 
V(»-i)2p+1  assoc'alcd  with  an  averaging  time  scale  of  2P+1  data  points,  and  the  sum  of 
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Dv(p+1)  is  the  total  signal  variance.  The  MR  spectra  have  units  of  variance  as  a 
function  of  averaging  time-scales  from  which  upper  and  lower  time-scale  limits  can  be 
obtained  and  converted  to  coarse  frequency  scales.  An  example  of  the  MR  response  to 
a  synthetic  periodic  time  series  signal  and  resulting  test  spectra  is  shown  in  Appendix 
C.  The  test  revealed  MR  limitations  that  affect  the  resolution  of  the  frequency 
boundaries,  such  that  MR  only  provides  coarse  limits  for  the  frequencies  of  the 
motions. 

The  averaged  MR  spectra  for  each  zone  are  shown  in  Figure  13.  The  32  Hz 
sampling  rate  of  the  ADV  allows  defining  the  higher  frequency  limit  of  the  MR 
analysis.  From  the  test  spectra  in  Appendix  C,  the  peak  in  the  MR  spectrum 
corresponds  to  the  frequency  that  captures  the  higher  frequency  limit  of  the  signal 
variance  or  energy.  The  MR  spectra  of  the  ADV  stream-wise  velocities  consistently 
have  peaks  at  0.125  Hz.  To  the  left  of  the  MR  spectra  peaks,  the  variance  steadily 
decreases  and  where  the  spectra  levels  out  to  the  lowest  frequency,  or  approaches  zero, 
indicates  the  lower  frequency  limit.  This  limit  consistently  occurs  at  0.001  Hz.  The 
periods,  T,  of  the  coherent  flow  structures  were  computed  using  Eq.  2  where  //  is 
0.0039  Hz  and/?  is  determined  from  energy  spectra.  Table  4  lists  the  variables  T,  L, 
and  L/h  for  each  dataset,  where  L  is  the  coherent  flow  structure  length  determined  from 
Eq.  1.  The  length  scales  range  from  25  to  91  m,  however,  care  must  be  exercised  when 
interpreting  the  scales  of  motion  detennined  from  this  method.  The  low-frequency 
limit  is  an  in  extremis  boundary,  that  is,  the  signal  contains  no  motion  less  than  0.001 
Hz,  but  the  actual  limit  is  between  0.001  and  0.0039  Hz.  Choosing  an  incorrect 
frequency  limit  will  skew  the  mean  T,  and  therefore,  the  lengths  in  Table  3  represent 
another  coarse  resolution  description  of  the  possible  lengths  for  coherent  flow 
structures. 
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IV.  DISCUSSION 


The  results  provide  an  indication  of  the  potential  length  scales  of  large-scale 
coherent  flow  structures  in  a  deeper  and  faster  gravel-bed  river  than  previously 
measured  in  the  field.  Energetic,  low-frequency  motions  containing  40  to  50%  of  the 
spectral  energy  were  found  for  all  deployment  locations.  The  low-frequency  motions 
are  identified  as  large-scale  coherent  flow  structures  as  measured  by  their  relatively 
long  coherent  length-scale.  Lc  for  the  low-frequency  structures,  which  averaged  25  m 
in  the  slower  flows  and  10  m  in  the  faster  flows  (Table  2),  were  shorter  than  their 
corresponding  motion  lengths  (Table  3).  This  implies  the  flow  structures  are  not 
coherent  for  their  entire  length,  and  suggests  the  structures  evolve  as  they  propagate 
downstream  that  is  shown  in  Figure  10  as  a  decrease  in  coherency  with  lag. 

Decay  in  coherence  is  caused  by  either  a  change  in  amplitude  or  a  shift  in  period 
as  the  structure  passes  through  the  array  (Bendant  and  Piersol  2000),  and  appears  to 
occur  more  rapidly  in  the  faster,  more  turbulent  flows.  There  are  several  possibilities 
responsible  for  this  process.  A  change  in  period,  or  frequency  shift,  may  be  the  result 
of  an  energy  transfer  between  scales  of  motions  as  noted  by  Carrasco  and  Vionnet 
(2004)  in  a  turbulent  flow  experiment.  If  the  coherent  flow  structures  are  considered 
alternating  regions  of  high  and  low-speed  wedges  as  described  by  Roy  et  al.  (2004),  Lc 
may  correspond  to  a  single  wedge.  The  succeeding  wedge  would  then  appear  as  a 
structure  with  different  frequency  or  amplitude. 

Dinehart  (1999)  discussed  the  potential  of  a  group  of  high-speed  pulses 
appearing  as  low-frequency  motions,  which  can  result  from  Fourier  analysis  creating  an 
apparent  period.  For  example,  IMFE  predicted  wavenumbers  are  based  on  fitting 
periodic  functions  to  the  velocity  co-spectra,  and  therefore,  the  method  “forces”  the 
results  to  be  periodic.  The  results  from  f-ks  indicate  that  the  ADCP  beam  length  was 
not  long  enough  to  reliably  capture  either  the  motion  lengths  or  the  coherent  lengths  of 
the  low-frequency  flow  structures.  A  longer  (>  20  m)  beam  length  is  needed  to  capture 
these  flows.  This  substantial  sized  array  explains  why  most  studies  have  been  confined 
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to  slow,  shallow  flows  requiring  shorter  arrays  and  questions  the  validity  of  a  2  to  6 h 
length  scale  for  fast,  deep  flows.  A  longer  array  will  only  confirm  the  linearly 
extrapolated  values  for  Lc,  but  longer  motion  lengths  may  still  remain  unresolved. 
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V.  SUMMARY  AND  CONCLUSIONS 


This  study  was  motivated  by  Dinehart’s  (1999)  observation  of  large  coherent 
flow  structures  in  a  relatively  fast  and  deep  gravel-bed  river  that  did  not  fall  within  a  2 
to  6 h  length  scale.  Several  days  of  stream-wise  river  velocity  measurements  were 
examined  in  a  braided  reach  section  of  the  Kootenai  River,  ID  with  the  objective  of 
detennining  a  characteristic  length  scale  for  large  coherent  flow  structures.  A  custom 
ADCP  head  was  validated  for  use  as  a  reliable,  accurate,  and  convenient  instrument  for 
obtaining  stream-wise  velocities  at  relatively  long  beam  lengths  with  finer  spatial 
resolution  and  less  influence  on  the  flow  than  a  comparable  ECM  array.  In  the  outer 
region  of  river  flow,  low-frequency  coherent  structures,  or  rolling  vortices  as  described 
by  Shvidchenko  and  Pender  (2001),  propagate  downstream  at  the  advective  velocity  of 
the  flow  according  to  Taylor’s  hypothesis  as  confirmed  by  f-ks  spectra.  The  structures 
have  coherent  lengths  from  5  to  22  m  (Table  2)  and  motion  lengths  that  could  be  as 
long  as  25  to  90  m  (Table  4).  Dinehart  (1999)  noted  similar  large  scales  that  were 
attributed  to  gravel-bed  form  migration  during  storm  flow  events,  neither  of  which 
occurred  in  the  flows  of  this  study.  Coherency  results  suggest  the  structures  evolve  in 
either  amplitude  or  period  prior  to  reaching  their  full-length  scale.  The  specific  change 
of  the  structures  was  not  detennined,  but  indications  are  that  some  mechanism  caused 
the  motions  to  alter  their  characteristics.  Although  this  study  did  not  provide  a  precise 
stream-wise  length  scale,  it  succeeded  in  showing,  through  spectral  techniques  unique 
to  river  applications,  that  these  structures  are  relatively  large,  have  long  lengths,  and  are 
coherent  for  long  distances.  The  study  also  shows  that  a  record  length  longer  than  5 
hours  is  required  to  resolve  the  low-frequency  coherent  motions  with  statistical 
confidence  and  an  array  length  of  20  m  or  greater  is  needed  to  capture  the  coherent 
lengths  of  large-scale  flow  structures  within  the  measured  flows.  In  summary, 
energetic,  low-frequency  motions  were  measured  in  fast  and  deep  flows.  These 
motions  are  coherent  along  the  length  of  the  array  with  an  Lc  determined  to  be  longer 
than  the  array.  The  coherent  flow  structures  are  longer  than  a  2  to  6 h  scale,  and  the 
structures  persist  shorter  than  their  motion  lengths.  These  observations  indicate  an  area 
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of  future  research  to  determine  the  life  cycle  of  coherent  flow  structures  to  increase 
understanding  of  their  generation  and  the  mechanism  that  alters  their  lengths  in  a 
relatively  short  distance. 
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APPENDIX  A:  CUSTOM  ACOUSTIC  DOPPLER  CURRENT 
PROFILER  VALIDATION 


An  increasing  number  of  riverine  environment  studies  have  used  ADCPs  to 
acquire  velocity  and  depth  measurements  (Gunawan  et  al.  2010;  Muste  et  al.  2004). 
These  instruments  are  relatively  inexpensive  and  easy  to  use.  A  typical  ADCP  uses 
three  angled  acoustic  beams  to  measure  velocities  which  pass  through  a  measurement 
volume  at  some  fixed  distance  from  the  instrument  head  and  requires  a  rotation 
calculation  to  convert  the  measurements  into  the  local  coordinate  axes.  The  2  MHz 
ADCP  used  herein  was  fitted  with  a  custom  head  that  uses  one  horizontal,  1 .7  degree 
beam  width  to  measure  along  beam  velocities  in  the  stream- wise  direction  at  equal  35 
cm  bins,  up  to  10  m  in  length,  with  a  sampling  frequency  of  1  Hz.  The  stream- wise 
oriented  beam  conveniently  serves  as  a  high-resolution,  equal-spaced  instrument  array 
that  is  non-intrusive  and  in  a  small  package,  however,  ADCPs  are  inherently  noisy 
instruments.  Also,  the  acoustic  beam  is  susceptible  to  side-lobe  scattering  by  surface 
and  bottom  interaction  if  the  narrow  ADCP  beam  width  is  not  properly  aligned  in  the 
horizontal.  The  velocity  measurements  collected  by  the  rigidly  mounted,  ADCP  head 
are  compared  to  the  ECM  array  and  ADV.  Since  the  ADCP  had  the  lowest  sampling 
frequency,  variances  were  computed  from  0.0039  Hz  to  the  Nyquist  frequency  of  0.5 
Hz.  The  difference  in  spectral  variances  between  the  ADCP,  ADV,  and  an  ECM  sensor 
was  consistently  less  than  10%  (e.g.,  Figure  14).  The  spectra  of  the  three  instruments 
were  equal  in  amplitude  as  a  function  of  frequency  and  are  not  aliased  by  the  coarse  1 
Hz  sampling  frequency.  The  f-ks  spectra  from  the  ADCP  and  ECM  array  also 
compared  well  with  similar  energy  and  wavenumber  estimates  (Figure  14).  The  results 
validate  the  ADCP  velocity  measurements  for  data  analysis  and  show  that  the  custom 
ADCP  is  a  reliable  and  convenient  instrument  for  resolving  the  scales  associated  with 
large  coherent  flow  structures. 
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APPENDIX  B:  ENERGY  SPECTRA  ESTIMATES 


Fourier  spectral  analysis  is  explained  in  detail  in  most  time  series  analysis  texts 
such  as  Bendant  and  Piersol  (2000)  and  Emery  and  Thomson  (2001).  To  compute 
energy  spectra,  averaging  is  done  over  a  number  of  windows  with  a  choice  of  window 
type.  A  Hanning  tapered  window  was  chosen  in  this  study  to  reduce  spectral  leakage 
and  it  provides  a  balance  of  resolution  and  statistical  confidence  (Biltoft  and  Pardyjak 
2009).  The  relationship  between  frequency  resolution  and  window  length  is  defined  by 

V  =  (B-D 


where  Tw  is  the  window  length.  A  window  length  of  256  s  (A/ =0.0039  Hz)  was 
chosen  as  a  compromise  between  frequency  resolution  and  total  record  length  and 
follows  Sukhodolov  and  Uijttewaal  (2010)  recommendation  for  a  window  length  on  the 
order  of  a  few  hundred  seconds  to  minimize  standard  error  and  variance  losses. 
Degrees  of  freedom,  which  is  the  number  of  statistically  independent  values  used  in 
each  spectral  estimate,  is  determined  for  a  Hanning  window  by 


DOF  =  — x 
3 


(B-2) 


where  Tt  is  the  total  record  length  (Emery  and  Thomson  2001).  Energy  spectra  of  each 
velocity  time  series  asymptote  to  a  common  106  DOF,  which  denotes  the  number  of 
windows  that  captures  most  of  the  signal  variance  (Lesht  1980).  The  selected 
parameters  result  in  7Yof  ~1.5  hours. 


A  significant  spectral  peak  was  found  in  the  middle  record  of  a  longer  ADCP 
dataset  (~  5  hour  record  length).  Despite  the  process  used  to  establish  ideal  Fourier 
analysis  parameters  just  described,  a  test  was  conducted  to  determine  if  a  smaller  Tt 
would  be  valid  for  analysis.  A  Zone  2  dataset  was  divided  into  Tt  sizes  of  10,  20,  30 
and  60  min  (Tw  =  256  s)  (Figure  15).  The  spectra  of  the  smaller  7Y  returned  significant 
peaks  suggesting  clearly  defined  periodic  motions  within  the  flow.  Even  though  DOF 
were  reduced  in  the  smaller  Tt  ,  the  peaks  are  significant  based  on  95%  confidence 
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intervals  computed  from  a  chi-square  test  value  (Emery  and  Thomson  2001).  A 
consistent  peak  across  consecutive  spectra  does  not  appear  and  the  frequency  band  at 
which  a  peak  occurs  is  not  consistent.  A  random  Gaussian  noise  signal  was  generated 
with  the  mean  and  standard  deviation  (1.4  and  0.23)  of  the  river  velocities  to  test  the 
validity  of  the  peaks  and  Tt  under  the  premise  that  any  signal  can  appear  periodic, 
given  a  small  enough  7>.  Surprisingly,  significant  peaks  are  noted  in  the  noise  signal 
spectra  (Figure  16)  up  to  a  ~1.5  hour  Tt.  This  result  suggests  that  any  significant  peaks 
in  energy  spectra  from  a  Tt  less  than  1.5  hours  cannot  be  confidently  discernable  from 
random  white  noise. 
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APPENDIX  C:  MULTI-RESOLUTION  DECOMPOSITION  TEST 

SPECTRA 


A  single  synthetic  periodic  signal  composed  of  several  low-frequency  periodic 
sub-signals  was  created  to  analyze  the  behavior  of  MR.  The  amplitudes  of  the  sub¬ 
signals  were  held  constant  since  no  appreciable  change  was  noted  in  the  MR  spectra 
when  the  amplitudes  were  varied.  The  range  of  frequencies  for  the  synthetic  signal  was 
chosen  to  produce  MR  spectra  similar  to  field  datasets.  A  spectrum  and  MR  of  the 
synthetic  signal  is  shown  in  Figure  17.  The  frequencies  of  the  synthetic  sub-signals  are 
plotted  on  the  MR  to  highlight  how  the  shape  of  the  resulting  MR  spectrum  is  used  to 
determine  signal  frequency  limits.  The  peak  in  the  MR  spectrum  corresponds  to  the 
averaging  time-scale,  which  captures  the  higher  frequency  limit  of  the  signal.  The  MR 
spectrum  decreases  from  the  peak  to  lower  frequencies  and  the  point  at  which  the 
spectrum  levels  out  or  approaches  zero  denotes  the  lower  frequency  limit  of  the  signal. 
The  test  revealed  MR  is  limited  by  a  significant  amount  of  spectral  leakage  associated 
with  the  2P  time-scale  separation,  and  the  variance  of  the  motions  for  a  particular  time- 
scale  are  shifted  to  the  next  lower  time-scale. 


27 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


28 


APPENDIX  D:  FIGURES  AND  TABLES 


Figure  1 .  Overview  of  a  section  of  the  Kootenai  River,  ID  and  instrument 
deployments  (circles).  Separation  of  deployments  into  zones  is 
highlighted  (boxed  areas)  and  described  in  Figure  2 
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Zone  1 


Figure  2.  Diagram  depicting  separation  of  zones  used  in  this  study.  Zone  1 
encompasses  a  section  of  the  main  channel.  Zone  2  is  located  2  km 
upstream.  Zone  3  includes  the  side  channel,  which  contains  a  riffle 
section.  Zone  4  is  located  downstream  of  a  man-made  concrete  piling  100 
m  upstream  of  Zone  3. 
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Figure  3.  Picture  of  frame  with  ECM  array  (not  visible),  ADV,  and  ADCP  deployed 
in  the  Kootenai  River,  ID. 


Figure  4.  ADCP  velocity  spectra  for  representative  locations  in  each  zone  (Figure 
2)  with  95%  chi-square  confidence  limits.  Separation  of  motions  between 
the  low-frequency  energetic  structures  and  turbulent  energy  cascade  is 
denoted  by  solid  vertical  lines. 
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Spectra  of  Zone  2  (Figure  2)  ADCP  velocities  computed  up  to  the  noise 
frequency  limit  (0.2  Hz)  using  varying  window  lengths  up  to  the  total 
record  length  (~1.5  hrs). 
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Figure  6.  Normalized  logarithmic  spectra  of  all  ADV  deployments  with  mean 
slopes  of  associated  motions  within  the  flow  (low-frequency  motions, 
mid-frequency  turbulence,  and  high-frequency  noise). 
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Normalized  Cumulative  Spectral  Density  [%] 


Figure  7.  Normalized  cumulative  spectral  density  for  representative  deployments  in 
each  zone  (Figure  2)  denoting  the  percentage  contribution  of  the  motions 
to  the  total  variance  per  frequency  bin,  disregarding  the  noise  floor  (f  > 
0.25  Hz). 
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Figure  8.  Coherency  determined  between  the  velocity  measurement  at  the  first  bin 
and  the  velocity  measurement  at  various  lags  up  to  the  maximum  ADCP 
beam  length  for  representative  deployments  for  Zone  2  (a)  and  Zone  4  (b) 
(Figure  2).  Solid  horizontal  line  denotes  significance  level  (0.06). 
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Figure  9.  The  bin  that  retains  the  maximum  coherent  distance  is  plotted  as  a 

function  of  frequency  for  representative  deployments  of  Zone  1  (solid 
line),  Zone  2  (dotted  line),  Zone  3(dashed  line),  and  Zone  4  (dashed-dot 
line)  (Figure  2).  Triangle  markers  indicate  coherence  exceeds  the 
maximum  ADCP  bin  distance. 


Figure  10.  Squared  coherency  function  computed  for  the  frequency  band  of  low- 
frequency  coherent  flow  structures  as  a  function  of  ADCP  bin  lag  for 
representative  deployments  in  each  zone  (Figure  2).  Coherent  lengths  are 
extrapolated  from  the  fitted  line. 
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Figure  11.  Normalized  co-spectra,  as  a  function  of  bin  lag,  at  various  frequencies  for 
ADCP  velocities  from  a  Zone  2  deployment  (Figure  2). 


Figure  12.  Frequency-wavenumber  spectra  of  velocity  time  series  (ADCP)  with 
associated  advective  velocity  (dotted  line)  according  to  Taylor’s 
hypothesis  (Eq.  1)  for  representative  deployments  in  each  zone  (Figure  2). 
Vertical  lines  represent  the  theoretical  maximum  motion  length  resolution 
of f-ks  as  either  two  (solid)  or  four  times  (dashed)  the  ADCP  beam  length. 


37 


Figure  13. 


Multi-resolution  decomposition  of  the  stream-wise  ADV  velocities.  MR 
spectra  for  each  deployment  are  averaged  for  each  of  the  zones  (Figure  2) 
with  standard  deviation  error  bars.  Vertical  lines  denote  frequency  limits 
of  the  signals. 
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Figure  14.  Spectra  (top)  of  ECM  (circles),  ADV  (triangle),  ADCP  (square)  velocity 
measurements  and/:/cv  spectra  (bottom)  for  a  Zone  2  deployment  (Figure 
2).  Spectral  variances  are  computed  up  to  the  ADCP  Nyquist  frequency 
(0.5  Hz).  Differences  of  variance  between  the  ADCP  (0.0179  m2/s2)  and 
the  ECM  (0.0180  m2/s2)  and  ADV  (0.0195  nr/s2)  were  computed  as 
0.78%  and  8.6%.  F-ks  spectra  (bottom)  of  ECM  and  ADCP  array 
compared  well. 
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Figure  15. 


ADCP  velocity  spectra  of  a  Zone  2  deployment  (Figure  2)  divided  into 
record  lengths  of  10  (a),  20  (b),  30  (c),  and  60  (d)  min  with  corresponding 
95%  chi-square  confidence  limits.  Significant  spectral  peaks  are 
highlighted  (circles). 


40 


Figure  16. 


Spectra  of  a  synthetic  random  Gaussian  noise  signal  computed  with 
different  record  lengths  of  10  (a),  30  (b),  60  (c),  and  90  (d)  min  with 
corresponding  95%  chi-square  confidence  limits.  Significant  spectral 
peaks  are  highlighted  (circles). 
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Figure  17.  Fourier  spectrum  (top)  and  MR  (bottom)  of  a  synthetic  signal  composed 
of  periodic  signals.  Vertical  dashed  lines  represent  corresponding 
frequencies  of  the  signals.  Solid  vertical  lines  indicate  frequency  limits 
determined  from  MR. 


Zone  #  /Yearday 
(2011) 

U  [m/s] 

h  [m] 

River  feature 

1/221 

1.52 

0.63 

Main  channel 

1/221 

1.38 

0.74 

Main  channel 

1/224 

1.63 

0.78 

Main  channel 

2/225 

1.28 

0.61 

Main  channel  -  upstream 

2/225 

1.47 

1.03 

Main  channel  -  upstream 

2/226 

1.41 

0.72 

Main  channel  -  upstream 

3/227 

0.59 

1.70 

Small  channel 

3/228 

0.46 

0.61 

Small  channel 

4/222 

0.57 

1.19 

Downstream  of  man-made  structure 

4/224 

0.29 

0.89 

Downstream  of  man-made  structure 

Table  1.  List  of  frame  deployments  and  corresponding  zones  (Figure  2)  used  for 
data  analysis  with  associated  mean  stream-wise  velocity,  U,  depth  of 
flow,  h,  and  major  river  feature. 
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Zone  # 

h  [m] 

Lc  [m] 

Tc  [s] 

L/h 

1 

0.63 

5 

3 

8 

0.74 

10 

7 

14 

0.78 

8 

5 

10 

2 

0.61 

8 

6 

13 

1.03 

7 

5 

7 

0.72 

10 

7 

14 

3 

1.70 

20 

33 

12 

0.60 

22 

48 

37 

4 

1.20 

7 

12 

6 

0.89 

22 

73 

25 

Table  2. 

Linearly  extrapolated  coherent  lengths,  times,  and  length  scales  as  a 
function  of  flow  depth,  h,  for  each  zone  (Figure  2). 

Zone  # 

ADCP  Beam  Length  [m] 

a)  L  [m]  (cs,  turb) 

b)  L  [m]  (cs,  turb) 

1 

*2.8 

? 

4.2 

? 

— ,  3.9 

4.9 

? 

2.4 

2 

5.25 

— ,  3.2 

4.7 

3.85 

? 

— ,  4.3 

6.3 

? 

— ,  6.5 

3 

*10 

4.7 

7.8,  4.7 

*9.5 

— ,  4.5 

6.5,  3.9 

4 

*2 

3.5 

2.1 

2.4 

Table  3.  Motion  lengths  computed  from  f-ks  spectra  using  two  times  the  ADCP 

beam  length  (a)  and  four  times  the  ADCP  beam  length  (b)  for  coherent 
flow  structures  (cs)  and  turbulence  (turb)  of  each  zone  (Figure  2). 
Asterisks  denote  the  beam  length  prevented  the  calculation  of  the  full 
range  of  wavenumbers  and  lines  indicate  the  wavenumbers  are  beyond 
resolvable  limits. 
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Zone  # 

h  [m] 

ns] 

L=UT  [m] 

L/h 

1 

0.63 

30 

45 

71 

0.74 

65 

91 

123 

0.78 

36 

58 

74 

2 

0.61 

35 

45 

74 

1.03 

54 

80 

78 

0.72 

40 

56 

78 

3 

1.70 

74 

44 

26 

0.60 

82 

38 

63 

4 

1.20 

84 

48 

40 

0.89 

87 

25 

28 

Table  4.  Length  and  temporal  scales  of  the  coherent  flow  structures  in  each  zone 
(Figure  2)  computed  from  Eqs.  1  and  2  using  frequency  limits  determined 
from  MR  and  energy  spectra. 
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